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ABSTRACT 

The  role  of  acoustics  in  flame/vortex  interactions  is  examined  via  asymptotic  analysis  and  numeri¬ 
cal  simulation.  The  model  consists  of  a  one-step,  irreversible  Arrhenius  reaction  between  initially 
unmixed  species  occupying  adjacent  half-planes  which  are  allowed  to  mix  and  react  by  convec¬ 
tion  and  diffusion  in  the  presence  of  an  acoustic  field  or  a  time-varying  pressure  field  of  small 
amplitude.  The  main  emphasis  is  on  the  influence  of  the  acoustics  on  the  ignition  time  and  flame 
structure  as  a  function  of  vortex  Reynolds  number  and  initial  temperature  differences  of  the  reac¬ 
tants. 
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1.  INTRODUCTION 


Acoustic  waves  play  a  pervasive  role  in  compressible  flows  even  in  the  low  Mach  number  limit.  They 
are  a  necessary  concomitant  of  nonlinear  interactions  such  as  those  present  in  turbulent  flows.  For  instance, 
the  interaction  of  any  disturbance,  whether  it  is  of  the  entropy  or  the  vorticity  type,  with  steep  gradients  as 
in  shocked  flows  [e.g.,  McKenzie  and  Westphal(I968),  Zang,  Hussaini  and  Bushnell  (1984)]  generates 
acoustic  waves.  Chemical  reactions  enhance  them  [Toong,  et  al.  (1974),  Roberts  (1978).  Clarke  (1985).  and 
Jackson,  Hussaini  and  Ribner  (1993)].  Conversely,  acoustic  waves  can  engender  instability  in  laminar  shear 
flows  which  may  eventually  suffer  transition  to  turbulence  [Goldstein  and  Hultgren  (1989)].  In  reactive 
flows,  they  play  a  similar  critical  role.  It  is  now  well  established  that  they  can  significantly  affect  the  stabil¬ 
ity  of  flames  [McIntosh  (1986)],  the  transition  to  detonation  [Urtiew  and  Oppenheim  (1966).  Lee  (1977), 
Shepherd  and  Lee  (1992)],  etc.  Oran  and  Gardner  (1985)  provide  a  fine  review  of  combustion-acoustics 
interactions  with  an  emphasis  on  the  physics  rather  than  the  mathematical  analysis  and  techniques. 

In  the  present  study,  the  role  of  acoustics  in  the  flame/vortex  interaction  process  is  investigated.  In 
order  to  provide  a  proper  perspective,  the  literature  on  the  flame/vortex  interaction  problem  is  briefly 
reviewed  here.  This  problem  was  first  formulated  by  Marble  (1985)  to  model  one  of  the  basic  physical 
mechanisms  underlying  the  complex  processes  of  turbulent  diffusion  flames  or  combustion  in  vortex- 
dominated  flows.  The  model  consists  of  a  diffusion  flame  with  fast  chemical  kinetics  (along  the  horizontal 
axis  separating  two  reactants  occupying  the  upper  and  lower  half-planes)  which  is  distorted  by  a  vortex  with 
its  center  at  the  origin.  The  theoretical  results  have  established  that  the  flame  sheet  is  rolled  up  into  a  spiral 
around  the  vortex  on  the  convection  time  scale  forming  a  reacted  core,  and  then  spreads  across  the  spirals 
on  the  diffusion  time  scale;  the  growth  rate  of  the  reacted  viscous  core  obeys  a  similarity  law  as  well  as  the 
reactant-consumption  rate  which  is  independent  of  time.  These  physical  features  are  confirmed  by  numeri¬ 
cal  simulations  [Laverdant  and  Candel  (1988),  Rehm,  Baum,  Lozier  and  Aronson  (1989).  and  Norton 
(1983)]  implying,  thereby,  that  the  inclusion  of  shear  (which  was  neglected  by  Marble)  or  finite-rate  chemis¬ 
try  does  not  alter  the  qualitative  picture.  The  premixed  combustion  in  a  vortex  also  displays  similar  features 
[Peters  and  Williams  (1988)]. 

The  aforementioned  studies  assume  the  simultaneous  existence  of  a  flame  and  a  vortex  and  the  subse¬ 
quent  evolution.  Thus,  they  preclude  the  ignition  regime  which  is  the  main  thrust  of  Maca*aeg.  Jackson  and 
Hussaini's  work  (1992).  Their  work  is  based  on  the  assumption  of  a  constant  density,  one-step  Arrhenius 
reaction  between  the  fuel  and  the  oxidizer,  occupying  adjacent  half  spaces,  in  the  presence  of  a  vortex  with 
its  center  on  the  axis  of  demarcation  between  the  reactants.  Their  study  focuses  on  the  ignition  time,  loca¬ 
tion  and  the  flame  structure  as  a  function  of  the  vortex  Reynolds  number  and  the  initial  temperature 
differences  of  the  reactants.  In  the  absence  of  the  vortex,  the  problem,  of  course,  reduces  to  the  classical 
case  of  Linan  and  Crespo  (1976)  who  analyzed  the  continuous  temporal  evolution  of  such  a  configuration 
from  nutrly  frozen  flow  to  near  equilibrium  flow.  The  reactants  mix  by  diffusion  until,  at  some  finite  time,  a 
thermal  explosion  occurs  at  a  well  defined  location,  and  ignition  takes  place.  After  ignition,  a  pair  of 
deflagration  waves  (or  premixed  flamelets)  emerge  according  to  classical  thermal  explosion  theory.  These 
premixed  flamelets  are  quite  weak  in  that  the  temperature  rise  associated  with  them  is  small,  and  they  exist 
only  until  all  of  the  deficient  reactant  is  consumed.  Just  beyond  the  deflagration  waves,  a  diffusion  flame 
regime  exists  where  the  mixing  process  is  governed  by  diffusion  in  the  direction  normal  to  the  flame.  As 
time  increases,  the.  diffusion  flame  approaches  a  flame  sheet.  We  note  that  the  existence  of  a  well  defined 
ignition  point  and  the  premixed  flamelets  depends  critically  on  the  relative  magnitudes  of  the  two  initial 
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temperatures  to  that  of  the  adiabatic  diffusion  flame  temperature.  That  is.  if  the  adiabatic  diffusion  flame 
temperature  is  greater  than  either  initial  temperature,  a  well  defined  ignition  point  always  occurs,  followed 
by  the  premixed  flamelets.  On  the  other  hand,  if  the  adiabatic  diffusion  flame  temperature  is  between  the 
two  initial  temperatures,  there  is  no  well  defined  ignition  point,  and  a  single  premixed  flame  merges 
smoothly  into  the  diffusion  flame.  This  configuration  of  a  triple-ll.ame  (or  tribrachial  flame)  appears  to  he  of 
such  a  fundamental  nature  that  it  also  exists  in  steady  supersonic  reacting  laminar  mixing  layers  [e.g..  Jack- 
son  and  Hussaini  (1988),  Grosch  and  Jackson  (1991)1  and  flames  propagating  into  a  nonunifonn  mixture 
le.g.,  Buckmaster  and  Matalon  (19X8);  Dold  (19X9.  1991)]. 

The  initial  presence  of  a  vortex  alters  the  picture  considerably  (Macaraeg,  Jackson  and  Hussaini 
(1992)].  First,  a  hot  spot  develops  within  the  viscous  core  of  the  vortex  and  evolves  into  an  almost  circular 
flame  which  grows  with  time  according  to  a  similarity  rule  and  so  does  the  reactant-consumption  rate,  as  in 
the  case  of  Marble.  Next,  the  Linan  and  Crespo  scenario  or  the  tribrachial  flame  configuration  must  evolve 
at  infinite  distance  to  the  right  and  left  of  the  origin.  The  diffusion  flames  emanating  from  this  configuration 
then  moves  toward  and  finally  merges  with  the  flame  in  the  reacted  core  region.  The  above  picture  is  only 
valid  for  the  case  when  the  adiabatic  flame  temperature  is  greater  than  either  of  the  two  initial  temperatures. 
When  the  adiabatic  flame  temperature  lies  between  the  initial  temperatures,  a  well  defined  ignition  point 
does  not  occur,  but  rather  the  premixed  region  merges  smoothly  into  the  diffusion  flame  region,  independent 
of  the  vortex  Reynolds  number.  This  scenario  is  consistent  with  the  Linan  and  Crespo  analysis. 

The  presence  of  an  acoustic  field  or  time-dependent  pressure  perturbation  affects  both  the  ignition  and 
the  flame  structure.  To  investigate  the  ignition  regime,  an  asymptotic  theory  is  developed  here  which  holds 
strictly  for  near-equal  temperatures  of  the  reactants.  For  low  frequencies  of  a  time  periodic  pressure  pertur¬ 
bation,  ignition  can  be  accelerated  or  significantly  delayed  depending  on  the  phase.  For  moderate  to  high 
frequencies  ignition  is  always  enhanced.  An  approximation  based  on  homogeneous  ignition  theory  is  given 
which  predict  these  trends  very  well,  and  therefore  is  proposed  as  a  good  engineering  approximation  for 
obtaining  the  ignition  times  for  the  more  complicated  flow  fields.  For  the  complete  time  evolution  from 
nearly  frozen  flow  to  near  equilibrium  flow,  an  asymptotic  analysis  based  on  small  heat  release  is  employed. 
In  this  limit,  direct  numerical  simulations  show  that  acoustics  have  little  effect  on  the  flow  field  in  the  vor¬ 
tex  core,  while  the  flow  outside  the  core  oscillates.  It  is  believed  that  viscosity  in  the  reacted  core  damps  out 
the  acoustic  held. 

In  the  next  section  the  problem  is  formulated.  Section  3  presents  an  analysis  of  the  ignition  regime, 
valid  only  for  near  equal  initial  temperatures  of  the  reactants,  using  a  combination  of  large  activation  energy 
(and,  hence,  large  Zeldovich  number)  asymptotics  and  numerics  for  several  different  pressure  profiles.  The 
pressure  profiles  chosen  are  either  a  linear,  pulse,  or  a  sinusoidal  function.  This  last  case  is  stressed  since  it 
represents  a  single  Fourier  component  of  a  more  general  pressure  disturbance.  Also  in  this  section,  the 
approximation  based  on  homogeneous  ignition  theory  is  presented.  Section  4  presents  selected  numerical 
results  of  the  continuous  evolution  from  nearly  frozen  flow  to  near  equilibrium  flow  for  the  sinusoidal  pres¬ 
sure  profile  under  the  assumption  of  small  heat  release.  A  description  of  the  diffusion  flame  regime  is  given 
in  Section  5  and  our  conclusions  are  presented  in  Section  6. 

2.  PROBLEM  FORMULATION 

In  this  section  the  problem  for  the  time  evolution  of  initially  unmixed  species  occupying  adjacent 
half-planes  which  .are  then  allowed  to  mix  and  react  in  the  presence  of  a  vortex  subject  to  pressure 
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disturbances  is  formulated.  The  nondimensional  equations  governing  this  field,  assuming  constant  viscous 
and  thermodynamic  properties,  are  given  by  Buckmaster  and  Ludford  (1982)  and  Williams  (1985), 


pT  =  l  +  8  Pc(t)  +  y  M2P(x,y,zj),  (2.1a) 

Pt  +  (P«).t  +  (pv)y  +  (pn')z  =  0,  (2.1b) 

p  [  u,  +  u  ux  +vu,  +  w  u,  ]  +  Px  =  Sc  ( V2u  +  y  ATX),  (2.1c) 

p  [  v,  +  u  vx  +  v  vy  +  w  Vj  ]  +  /%  =  Sc  ( V-v  +  y  Ky ),  (2.  Id) 

p  f  w,  +  1 i  wx  +  v  wy  +  w  wz]  +  Pz  =  Sc  ( V2w  +  y  Kz ),  (2.  le) 


p[r,  +  „r.  +  vr,  +  *r.]-8^!^- 

-(Y-  1  )M1[P,  +  uPx  +vPy  +  w  PZ]  =  V2T  +  (y-  l)M2Sc  O  +  PS2.  (2.10 

p[F;.,  +uFJtX  +  vFj.y  +  wFjt2]  =  V2Fj  -ii,  j  =  1,2,  (2.1g) 

il  =  DapFiF2e'z>IT ,  (2.1h) 

where  K  =  ux  +  vy  +  w2 ,  O  is  the  viscous  dissipation  term,  and  V2  Is  the  three-dimensional  Laplacian 
operator.  Here,  (u,v,w)  are  the  velocity  components  in  the  (x,y ,z)  directions,  respectively;  p  is  the  den¬ 
sity;  T  is  the  temperature;  and  F  *  and  F2  the  mass  fractions  of  the  fuel  and  oxidizer,  respectively.  The  total 
pressure  P  =  1  +  8 Pc  +  yMlP  is  written  as  the  sum  of  the  compressible  component  Pc  plus  incompressi¬ 
ble  component  P ,  as  is  consistent  with  the  small  Mach  number  approximation  [e.g.,  Majda  (1984)],  with  8 
the  amplitude  of  the  compressible  component.  The  actual  size  of  8  will  be  chosen  in  the  course  of  the 
analysis.  The  chemical  model  is  assumed  to  be  a  one-step,  irreversible  Arrhenius  reaction.  The  nondimen¬ 
sional  parameters  appearing  above  are  the  Schmidt  number  Sc  =  v  /  D  assumed  equal  for  both  species  with 
D  the  species  diffusion  coefficient  and  v  the  kinematic  viscosity;  the  Zeidovich  number  Ze  -E  t(R°  T„) 
with  E  the  dimensional  activation  energy  and  R°  the  universal  gas  constant;  the  Damkohler  number  Da 
defined  as  the  ratio  of  the  characteristic  time  scale  id  to  the  characteristic  reaction  time  scale  iR ;  y  the  ratio 
of  specific  heats;  and  finally  P  the  heat  release  per  unit  mass  of  .  In  specifying  the  Damkohler  number 
in  the  course  of  the  analysis,  the  appropriate  choice  for  the  characteristic  time  scale  is  made.  The  density, 
temperature,  pressure  and  mass  fractions  were  nondimensionalized  by  their  initial  values  p„.  /„, 
P „  =  p mR°  T„  and  F1  „  ,  respectively,  that  would  prevail  if  there  were  no  pressure  disturbances.  The 
Lewis  number  was  assumed  to  Lc  unity,  which  implies  that  Sc  =  Pr .  where  Pr  =  p„  Cp  v  /  X  is  the  Prandtl 
number,  X  is  the  thermal  conductivity  and  Cp  is  the  specific  heat  at  constant  pressure.  Lengths  and  veloci¬ 
ties  are  referred  to  the  relevant  diffusion  characteristic  scales  ld  -  D  id  and  Ud  =  ldl  td ,  respectively.  In 
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this  case,  the  Reynolds  number  based  on  this  choice  of  / d  and  Ud  is  given  by  Re  =  ld  Ud  /  v  =  1/Sc.  The 
Mach  number  M  is  defined  as  Ud  / a„,  the  ratio  of  the  characteristic  diffusion  speed  to  the  speed  of  sound. 
The  above  choices  for  td,  ld,  and  Ud  are  consistent  with  the  scales  chosen  in  the  absence  of  a  vortex. 


An  analysis  of  the  ignition  zone  with  acoustic  interactions  must  take  into  account  the  two  fundamental 

ratios 


characteristic  time .  td  ld  /  Ud 

acoustic  time ,  ta  lalaM 


la_ 
Id ' 


These  two  ratios  are  related  via 


1 

N  M  ’ 


where  M  is  the  Mach  number  defined  above.  Here,  our  analysis  is  based  on  the  assumption  that  t  =  O  ( 1 ), 
and  hence  the  acoustic  wavelength  is  much  longer  than  the  characteristic  diffusion  length  scale  for  small 
Mach  numbers.  Thus,  the  pressure  gradient  within  the  combustion  zone  is  negligibly  small  and  the  pressure 
is  essentially  a  function  of  time,  and  so  the  flame  is  treated  as  isobaric  to  an  excellent  approximation.  In 
addition,  our  analysis  further  assumes  that  the  magnitude  of  the  acoustic  time  relative  to  the  reaction  time  is 
given  by 

characteristic  time ,  td  ta 

Da  = - - -  st  —  . 

reaction  time ,  tR  tR 


Taking  1  -  0(1),  we  shall  assume  that  the  reaction  time  is  much  faster  than  the  acoustic  time,  and  so  we 
take  the  Damkohler  number  to  be  proportional  to 


where  Ze  >>  1  is  the  Zeldovich  number  defined  above,  thus  defining  the  characteristic  time.  The  exact  form 
is  chosen  in  the  course  of  the  analysis.  In  our  analysis,  if  the  Mach  number  M  was  not  chosen  small,  we 
would  also  need  to  consider  spatial  variations  of  the  pressure  disturbance  inside  the  reaction  zone. 

McIntosh  (1986,  1989,  199Ia,b,  and  the  references  cited  therein)  pointed  out  the  importance  of  the 
two  fundamental  ratios  given  above  when  the  characteristic  velocity  is  given  by  the  flame  speed  Uf.  For 
this  case,  three  distinct  cases  of  flame-acoustic  interactions  arise  depending  on  the  magnitude  of  T.  and  these 
have  been  classified  in  McIntosh  (1991a).  Further  phenomena  can  exist  depending  on  the  magnitude  of  the 
acoustic  time/reaction  time  ratio,  and  a  discussion  of  this  can  be  found  in  Clarke  (1985).  For  x-  0(1), 
McIntosh  and  Wilce  (1991b)  point  out  that  audible  sound  and  its  interference  with  flames  would  come  under 
this  category,  and  a  considerable  body  of  work  exists  using  this  assumption  (e.g.,  references  cited  therein). 

We  shall  now  assume  that  the  flow  is  independent  of  the  downstream  direction  x  for  all  time. 
Although  not  true  for  any  real  flow  situation,  it  is  consistent  with  experimental  observations  that  streamwise 
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vortical  structures  in  turbulent  mixing  layers  are  essentially  aligned  with  the  flow  and  are  of  great  extent. 
Also,  neglecting  the  x  -dependency  can  be  thought  of  as  a  local  approximation  for  the  vortices  found  in  tur¬ 
bulent  flames.  This  type  of  approximation  has  also  been  used  previously  in  describing  the  evolution  of  a 
shear  flow  with  an  embedded  streamwise  vortex  [Corcos  (1988),  Pearson  and  Abernathy  (1984)}.  We  shall 
also  assume  that  the  Mach  number  is  small,  as  is  typical  for  flames  [e.g.,  Buckmaster  and  Ludford  (1982), 
Williams  (1985)].  Thus,  the  governing  equations  (2.1),  in  the  limit  of  small  Mach  number  and  neglecting 
all  dependency  in  x ,  reduce  to 


pr  -  l  +8PAD. 

(2.2a) 

P;  +  (pv)y  +  (pH-),  =  0, 

(2.2b) 

p  ( u,  +  v  Uy  +  w  uz  ]  =  Sc  V2« 

(2.2c) 

p[v,  +  v  v,,  +  w  vz  ]  +  Py  =  Sc  (VH'  +  J Ky), 

(2.2d) 

p [ w,  +  v  wy  +  w  wz]  +  P.  =  Sc  ( V2w  + 

(2.2e) 

pr r,  +  v7,  — =  v2r  +  pn, 

y  at 

(2.20 

p  [Fjj  +  v  Fj  y  +  w  Fj_z  ]  =  V2F;  -  £2,  j  =  1.2, 

(2.2g) 

where  V2  =  (  )yy  +  ( )zz  is  the  two-dimensional  Laplacian  operator,  K  =  vy  +  wz,  and  we  have  considered 
the  distinguished  limit  M 2  <<  8  <<  1.  Since  we  are  interested  in  reacting  flows  with  vortical  motion,  we 
transform  (2.2)  into  cylindrical  coordinates  (y ,  z )  -» (r ,  0)  by  using  the  transformations  y  =  r  cos  0  and 
z=rsin0.  If  we  define  g  to  be  the  radial  velocity  and  h  to  be  the  tangential  selocity,  i.e., 
v  =  g  cos  8  -  h  sin0  and  w  =  g  sin0  +  h  cos0,  then  (2.2)  transforms  to 

pr  =  1  +8Pc(r), 

(2.3a) 

P(  +  (  P  £  )r  +  ~(p^)e+  ~  0. 

(2.3b) 

p  [  U,  +  g  Ur  +  —  ue  ]  =  Sc  V2« , 

(2.3c) 

PU,  +  g  gr  +  j  g«  ~  -y- 1  +  Pr  =  Sc 


v2*  -4 


4/,«+ 

3 


(2.3d) 


-6- 


P  {h,  +  g  hr  +  ~  h0  + 

r  r 


W2h  - 


2 

~J  So  + 


(2-3e) 


p[T,  +  g  Tr  +  —Tb]  -8^ — =  V2T  +  pi2,  (2.30 

r  y  at 

P  [Fj.,  +  g  Fj  r  +  y  Fy.e]  =  V2Fj  -  f2,  /  =  1,2,  (2.3g) 

where  V2  is  now  the  two-dimensional  Laplacian  operator  in  cylindrical  coordinates,  ;tnd 
K  =  gr  +  g  / r  +  hBlr.  The  system  (2.3)  must  be  solved  subject  to  appropriate  boundary  and  initial  condi¬ 
tions. 


Since  it  is  known  that  small  pressure  disturbances  of  0(Ze  '*)  can  cause  0(1)  changes  in  certain  flame 
properties  such  as  the  burning  rate,  extinction,  transition  from  laminar  to  turbulence,  etc.  (Buckmaster  and 
Ludford  (1982),  Buckmaster  (1992),  Kapila  (19921],  we  shall  limit  our  attention  to  small  pressure  distur¬ 
bances  by  defining  for  the  total  pressure 

P  =  1  +  Ze"1  -3—  Pc  (t )  +  O  (A/2),  8  =  — 1—  Ze~l\  (2.4) 

Y-  1  Y"  1 

the  factor  y/(y-  1)  in  the  definition  of  8  was  chosen  for  convenience.  Thus,  we  are  concerned  with 
0(Ze~])  acoustic  perturbations  superimposed  on  a  basic  state.  The  pressure  perturbation  Pr(t)  can  be  deter¬ 
mined  by  solving  the  acoustic  equations  in  the  outer  field  fe.g..  Van  Harten,  Kapila  and  Matkowsky  (1984)], 
and  indicates  that  the  pressure  can  be  controlled,  for  example,  from  the  surroundings.  A  general  expression 
for  Pc  can  be  found  in  Majda  (1984)  and  will  depend  on  the  type  of  conditions  imposed  on  the  geometry 
under  consideration.  Since  we  are  interested  in  the  influence  of  pressure  variations  on  the  flame,  we  will 
not  solve  this  outer  problem  and,  henceforth,  consider  the  pressure  disturbance  Pr  to  be  prescribed. 

We  note  here  that  it  is  the  density  which  couples  the  momentum  equations  to  the  temperature  and 
mass  fraction  equations.  As  such,  the  numerical  solution  to  the  full  system  (2.3)  is  a  formidable  task.  One 
approximation  which  will  simplify  (2.3)  to  a  more  tractable  system  is  the  assumption  of  small  heat  release 
[e.g..  Van  Harten,  Kapila  and  Matkowsky  (1984)].  The  resulting  system  can  then  be  solved  numerically,  and 
will  be  presented  in  Section  4.  For  the  ignition  problem,  the  assumption  of  small  heat  release  is  not  needed 
in  order  to  simplify  the  system,  provided  we  restrict  the  analysis  to  near  equal  initial  temperatures,  and 
hence,  densities,  of  the  two  gases.  This  analysis  is  presented  in  the  following  section. 

3.  IGNITION 

At  time  t  =  0,  the  reaction  rate  is  exactly  zero  owing  to  the  product  F  ( F2  -  0.  For  /  >  0,  the  fuel  and 
oxidizer  begin  to  mix  by  diffusion,  as  well  as  by  convection  due  to  the  presence  of  the  vortex,  and  the  reac¬ 
tion  rate  is  no  longer  zero.  For  small  time,  it  can  thus  reasonably  be  assumed  that  the  effect  of  the  reaction 
on  the  overall  flow  field  is  small.  The  solution  for  small  time  in  which  the  reaction  rate  term  is  neglected  is 
known  as  the  inert  or  chemically  frozen  solution,  and  will  be  denoted  by  the  superscript  I.  In  what  follows, 
we  only  consider  the  case  of  near  equal  initial  temperatures.  This  assumption  leads  to  the  inert  solutions 
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T1  =  p/  =  1  +  0(Ze~‘),  where  the  0(Ze~')  term  is  included  to  allow  lor  small  initial  temperature 
differences  that  might  exist  in  the  absence  of  any  pressure  disturbances.  Thus,  the  leading-order  inert  solu¬ 
tion  corresponds  to  the  constant  density  approximation  prevalent  in  the  combustion  literature.  As  time 
increases,  more  of  the  combustible  mixes  until,  at  some  finite  time,  a  thermal  explosion  occurs  characterized 
by  significant  departure  from  the  inert.  To  analyze  the  ignition  process,  we  determine  the  effect  of  the  grow¬ 
ing  reaction  rate  by  expanding  about  the  inert  solution  as 

T  =  l  +  Ze~l[Ti-i-rcU)l  +  0(Ze^h  p  =  1  +  Ze~x  p>  +  0(/,e~2).  Fj  =  Fj  +  0(Ze~'),  (3.1) 


u  =u‘  +  0(Ze-‘),  g  =  g1  +  0(Ze').  h  =  h1  +  0(Ze "*), 


and  take  the  asymptotic  limit  Ze  — >  <*>.  The  leading-order  equations  are  given  by 


P.  +7-1  = 


r  r 


(3.2) 


u{  +  g1  uj  +  —  «e  =  Sc  W , 


g!  +  g1  g'r  +  -y  *e  -  ~  +  Fr  -  Sc 


V2g' 


r~ 


h{  +  g‘  ^  +  —  h's  +  sLhL  +  ipe  =  Sc 


W  -  —  +  — 


r 


Ti.,  +g'TUr  +  ~  T]  e  =  V2T,  +  A  F\  F'2eT)  +  PA‘\ 


(3.3a) 


(3.3b) 


(3.3c) 


(3.3d) 


(3.3e) 


(3.30 


Fj.,  +  g‘  Fj,r  + 


h' 


Fj,e  —  V2Fj,  j  =  1,2, 


(3.3g) 


where  we  have  chosen  the  Damkohler  number  Da  to  be 

Da a4) 

and  A  is  some  constant  that  will  be  chosen  in  the  course  of  the  analysis.  This  particular  choice  of  the  Dam¬ 
kohler  number  ensures  that  a  distinguished  limit  exists,  in  that  the  reaction  rate  term  is  of  the  same  order  in 
a  Zeldovich  number  expansion  as  the  time  derivative  terms  [see,  e.g.,  Buckmaster  and  Ludford  (1982)). 
Note  that  the  pressure  disturbance  now  appears  in  the  exponential  and  thus,  will  play  a  crucial  role  in  the 
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ignition  process. 

We  first  see  that  the  continuity  and  momentum  equations  (3.3b-e)  are  decoupled  from  the  rest  of  the 
system  and  thus,  will  completely  determine  u1 ,  g1 ,  h1  and  F.  At  t  -  0,  we  situate  a  point  vortex  at  the  ori¬ 
gin  and  allow  it  to  diffuse  under  the  action  of  viscosity.  Then,  assuming  the  flow  in  the  (r,0)  plane  to  be 
axisymmetric,  the  induced  flow  field  is  given  by  [Lamb,  1932] 

g'  =0,  h1  =  -  exp(~r2/4t  Sc)],  (3.5) 


P(r.t)  =  -f  —  dr.  (3.6) 

r 

and  u1  is  found  from  the  equation 

U I 

u\  +  —  u'B  =  Sc  V2m/.  (3.7) 

r 


The  solution  (3.5)-(3.6)  is  the  sotution  for  the  incompressible  Oseen  vortex,  where  R  =  T/2nv  is  the  vortex 
Reynolds  number  and  T  is  the  circulation  parameter.  The  vortex  Reynolds  number  R  is  considered  as  a  free 
parameter  and  its  choice  determines  the  importance  of  the  vortex  in  the  initial  mixing  of  the  reactants. 
When  R  =  0  (i.e.,  T  =  0),  there  is  no  vortex  and  the  initially  separated  reactants  mix  by  diffusion  only. 
When  the  vortex  Reynolds  number  R  is  non-zero,  the  initially  separated  reactants  mix  by  diffusion  and  con¬ 
vection.  The  solution  to  (3.7)  has  been  discussed  previously  in  the  context  of  streamwise-spanwise  vorticity 
interactions  for  nonreacting  flows  [Corcos  (1988),  Pearson  and  Abernathy  (1984)].  Note  that  «'  does  not 
appear  in  the  temperature  or  mass  fraction  equations  since  we  neglected  variations  in  x  coupled  with  the 
zero  Mach  number  approximation,  and  so  the  analysis  is  valid  for  either  streamwise  or  spanwise  vortices. 
The  velocity  component  ul  will  influence  the  temperature  equation  if  true  compressibility  effects  are  taken 
into  account;  such  will  necessarily  occur  after  ignition  has  taken  place. 

With  g1  and  h1  now  known,  the  following  equations  can  be  solved  for  T ,  and  F- 


i.< 


7-r,.  9  =  V' 


T  i  +  A  F\  F'j  e 


T,*Pe<f) 


(3.8) 


F'j.,  +  -7-^.0  =  V2/--;.  y  =  l,2.  (3.9) 

The  appropriate  boundary  and  initial  conditions  are  given  by 

Tj  =  0,  F\  =  1,  F\  =  0  at  t  =  0,  r  >  0,  0  <  8  <  it,  and  t  >  0,  r  -4  <»,  0  <  0  <  Tt, 

/r 

T\  -  Pr-  F[  =  0.  f'2  =  ~rr—  s  at  t=0,  r>0, 7t<0<27t, 

r  1  „ 


and  t  >  0,  r  »,  n  <  8  <  2  n. 


(3.10) 
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where  <}>  is  the  equivalence  ratio  defined  as  the  ratio  of  the  initial  mass  fraction  of  the  fuel  F  { to  the  initial 
mass  fraction  of  the  oxidizer  and  fl7  is  the  parameter  which  ttllows  for  small  initial  temperature 

differences  in  the  absence  of  pressure  disturbances.  If  0=  1,  the  mixture  is  said  to  be  stoichiometric;  if 
<t>  >  I  it  is  fuel  rich;  and  if  <j>  <  1.  it  is  fuel  lean.  Also,  if  j}r  is  less  than  zero,  the  oxidizer  (species  2)  is 
relatively  cold  compared  to  the  fuel  (species  1);  and  if  (37  is  greater  than  zero,  it  is  relatively  hot.  As  i 
increases,  the  solution  for  Tx  becomes  unbounded  at  some  finite  time  (t,K)  and  location  (y,f ,  r,? ).  This 
characterizes  the  ignition  regime.  The  special  case  R  -  0,  in  which  two  initially  unmixed  species  are 
allowed  to  diffuse  without  the  mixing  generated  by  the  vortex,  will  be  discussed  in  subsection  3.1,  while  the 
case  for  R  >0  will  be  presented  in  3.2.  Finally,  with  T j  determined,  the  density  perturbation  pi  can  be 
found  from  (3.3a). 

3.1.  R  -  0.  In  the  absence  of  any  vortical  motion,  the  governing  equations  reduce  to  the 
reactive -diffusion  equations 


Ti.,=Tu„ 

,  c{  r- /  ^  |  +  Pc(0 

+  A  F\  F  2  e 

(3.11) 

F1  - 
rj.i  ~ 

Fj,:z>  y=  1.2, 

(3.12) 

subject  to  the  boundary  and  initial  conditions 

T  i  =  0,  F\  -  1,  F\  =0 

at  l  =  0,  z  >  0,  and 

t  >  0,  z  — > 

(3.13a) 

=  0,  F\  =r‘ 

at  t  =  0,  2  <  0,  and  t 

>  0,  Z  — *  —  oo. 

(3.13b) 

The  parabolic  equations  for  the  mass  fractions  can  be  integrated  to  yield 

F\  =  T, 

* 

i 

7 

-©- 

II 

— .fN 

(3.14) 

where 


V  =  +  erf  t\), 

and  q  is  the  similarity  variable,  defined  as 


*1  = 


2  'll  ' 


(3.15) 


(3.16) 


Upon  substituting  the  inert  solutions  into  the  temperature  perturbation  equation,  and  transforming  to  the 
similarity  plate,  a  single  equation  for  T  j  emerges,  given  by 

4 1  T ,  (  -  2qT]  „  =  4t  T(1  -  V)eT'  (3.17) 

which  must  be  solved  numerically  subject  to  the  boundary  and  initial  conditions 
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T  |  =  0  at  i  -  0.  z  >  0.  and  T  >  0,  (3. 18a) 

7',  -  Pr  at  /=(),*<  0,  and  T  >  0,  z  -»  -■«>.  (3.1Kb) 

Here,  (he  constant  A  appearing  in  the  definition  of  the  D.imkohler  number  has  been  chosen  to  be  $  tor  con¬ 
venience.  The  special  case  of  no  pressure  disturbance  (Pc  =  0)  corresponds  to  the  results  of  Linan  ;tnd 

Crespo  ( 1 976)  which  has  been  discussed  in  the  Introduction.  Thus,  the  results  of  this  subsection  con  be 

thought  of  as  an  extension  of  the  Linan  and  Crespo  problem  fo  include  pressure  disturbances  In  the  follow¬ 
ing  subsections  we  will  select  three  different  pressure  variations,  and  compare  with  the  results  of  Lintin  and 
Crespo.  The  first  two  choices  for  the  pressure  variation  presented  below  corresponds  to  the  choices  taken  by 
Tedder  and  Kapila  (1991)  in  their  study  of  the  response  of  premixed  flames  to  pressure  perturbations.  The 
last  choice  will  be  a  sinusoidal  pressure  variation,  which  represents  a  single  Fourier  component  of  a  more 
general  pressure  disturbance. 

3.1.1.  UNEAR  PRESSURE  DECREASE  Let  the  flow  be  subject  to  a  linear  pressure  decrease  at 
t  =  0  of  the  form 

RrU)  =  -0)/.  (3.19) 

where  1  /  to  >  0  charactererizes  the  time  scale  of  the  pressure  disturbance.  The  influence  of  this  linear  pres¬ 
sure  profile  on  premixed  flames  has  been  discussed  by  Ledder  and  Kapila  (1991).  in  which  extinction  was 
seen  to  occur  for  all  values  of  to.  This  profile  w»  chosen  for  (his  study  because  of  its  simplicity  in  form  in 
that  the  ignition  lime  and  location  will  only  be  a  function  of  the  parameters  to  and  |37 .  and  also  may 
represent  the  first  term  in  a  Taylor  Series  expansion  of  a  more  general  pressure  disturbance.  Figure  1  is  a 
plot  of  the  maximum  of  T ,  versus  t  for  several  values  of  tu  and  p7  =  0.  The  case  to  =  0  corresponds  lo  that 
of  Linan  and  Crespo  (1976)  in  which  no  pressure  disturbance  is  present.  As  to  is  increased  from  zero,  the 
ignition  time  also  increases  until,  at  some  finite  value  to  s  tur  =  0.1395.  the  ignition  time  is  pushed  to 
infinity.  Thus,  small  values  of  u)  can  significantly  delay  ignition.  For  values  of  to  greater  than  tor .  the 
maximum  of  /',  initially  increases,  reaches  a  maximum,  then  begins  to  decrease  back  to  zero  as  time 
increases  further.  Thus,  ignition  does  not  occur  in  the  classical  sense.  Mathematically,  when  (he  pressure 
disturbance  becomes  large  and  negative,  the  nonlinear  source  term  of  (3.17)  becomes  exponentially  small 
and  so  can  be  neglected:  T ,  then  evolves  according  to  a  homogeneous  heat  equation.  Although  the  solution 
of  this  system  suggests  that  ignition  will  not  occur  for  to  grealer  loan  wr ,  we  note  that  the  presence  of  ;ui 
infinite  amount  of  fuel  and  oxidizer  necessitate  that  ignition  will  eventually  hike  place.  Asymptotically,  the 
expansion  (2.4)  breaks  down  when  t  is  of  the  order  0{Ze),  and  a  new  expansion  must  be  employed,  one 
which  allows  0(1)  changes  in  the  pressure  and  the  temperature;  the  waves  now  cease  to  be  acoustic  in 
nature.  The  leading-order  system  in  this  case  is  no  longer  that  with  constant  density,  and  hence,  the  full  sys¬ 
tem  (2.3)  must  be  solved  numerically.  This  numerical  problem  is  beyond  the  scope  of  the  present  study. 
Figure  2  is  a  plot  of  the  ignition  time  versus  to  for  three  values  of  the  temperature  ratio  parameter  pr. 
The  case  of  pT  =  0  was  presented  in  a  slightly  different  manner  in  Figure  1 .  Here,  it  is  clear  that  as  to 
approaches  0.1395  from  below,  the  ignition  time  goes  to  infinity,  and  for  to  greater  than  0.1395  there  is  no 
ignition.  The  case  m  =  0  corresponds  to  that  of  Linan  and  Crespo  (1974);  ignition  time  increases  as  p7  is 
decreased  from  zero  (/;2„  is  relatively  cold  compared  to  /•  |  „),  and  decreases  as  p7  is  increased  from  zero. 
For  each  value  of  the  temperature  ratio  parameter,  a  critical  value  of  to  is  seen  to  exist;  this  value  being  less 
than  0.1395  for  negative  pr  and  greater  than  0.1395  for  positive  values  of  pr.  In  the  absence  of  any 
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pressure  disturbance.  Linan  and  Crespo  have  shown  that  ignition  takes  place  along  a  line  parallel  to  the  v- 
axis  which  is  located  at  z,g  -  0  lor  =  0  and  resides  in  the  hotter  region  for  (i;  *  o.  In  the  presence  of  the 
linear  pressure  disturbance  ;tnd  for  Pr  -  0.  ignition  always  occurs  along  the  line  ru.  -  0  for  any  value  of  w. 
However,  for  pr  /  0,  the  ignition  line  z  ~  zlg  (c,s  <  0  for  |ir  >  1).  and  ztg  >  l)  for  P,  <  0)  moves  towtird  the 
line  z  =  0  as  to  is  increased. 

3.1.2.  PRESSURE  PULSE  Let  the  flow  be  subject  to  a  pressure  pulse  at  t  =  0  of  the  form 

PrU )  -  4Pa  ( 2~b”:  ~  ),  (3.20) 

where  aT1'2  >  0  is  the  pulse  width  and  |  PA  j  is  the  amplitude  of  the  pressure  pulse.  This  corresponds  to  a 
smooth  increase  beginning  at  zero,  increasing  or  decreasing  to  tin  extremum  \PA\  at  tor  =  1.  and  finally 
returning  to  zero.  The  influence  of  this  pressure  pulse  on  premixed  flames  has  been  discussed  by  Tedder 
and  Kapila  (1991),  in  which  extinction  is  seen  to  occur  for  both  negative  and  positive  amplitudes,  depending 
upon  the  parameters  of  the  problem.  Figure  3  is  a  plot  of  ignition  time  ilg  versus  to  for  vttrious  values  of  the 
amplitude  PA  with  pr  =  0.  The  case  of  PA  =0  corresponds  to  that  of  Linan  and  Crespo  (1976)  and  is 
shown  in  the  figure  as  a  dashed  line  for  reference.  In  addition,  all  data  curves  shown  in  this  figure  begin  at 
tig  -  5.816  at  0)  =  0,  again  consistent  with  the  results  of  Linan  and  Crespo.  For  fixed  PA  <  0,  the  ignition 
time  initially  increases  as  to  increases  from  zero,  reaches  a  maximum,  and  then  begins  to  decrease  as  to  is 
increased  further  reaching  an  asymptote  as  (0  ->  <».  For  fixed  to.  decreasing  the  pressure  amplitude  from 
zero  increases  the  ignition  time.  Thus,  ignition  is  delayed  for  negative  pressures.  In  contrast  to  the  previous 
case,  for  a  positive  pressure  pulse,  the  ignition  time  decreases  monotonically  as  to  is  increased  from  zero 
and  asymptotes  to  a  finite  non-zero  value  its  to  -»  hence,  ignition  is  enhanced.  Since  [V  =  ignition 
takes  place  along  a  line  parallel  to  the  y-axis  located  at  zig  =  0  for  any  vtilue  of  PA  and  to.  The  general 
behaviour  described  in  Figure  3  for  pr  =  0  is  seen  to  hold  for  tiny  value  of  [lr.  the  only  exception  being  the 

starting  value  of  tlg  at  to  -  0;  for  pr  <  0.  the  starting  vtilue  is  above  that  given  in  the  figure,  while  for 

Pr  >  0  it  is  below.  In  addition,  the  location  of  the  ignition  line  resides  in  the  hotter  region  and  moves 
toward  z  =  0  a>  to  is  increased  from  zero  for  fixed  PA . 

An  interesting  feature  is  seen  to  occur  as  PA  is  further  decreased  below  -  2.  This  is  shown  in  Figure  4 
where  we  plot  tlg  versus  to  for  PA  =  -4  and  Pr  =  C.  The  ignition  time  increases  slowly  up  to  about 
to  =  0.0018,  where  a  sudden  jump  in  tig  occurs.  As  to  increases  still  further,  the  ignition  time  reaches  a 

maximum  and  then  begins  to  decay  slowly,  reaching  an  asymptote  as  to  -»  ».  To  better  explain  this  jump 

phenomenon,  we  plot  in  Figure  5  the  maximum  of  T i  versus  time  for  severtil  values  of  to.  Note  that  there 
is  a  sudden  jump  in  the  ignition  time  between  to  =  0.0018  and  to  =  0.002.  For  vtdues  of  to  below  0.0018, 
the  maximum  temperature  perturbation  is  seen  to  increase  monotonically.  becoming  infinite  as  t  approaches 
tig.  For  values  of  to  greater  than  about  0.(X)18,  the  maximum  temperature  perturbation  increases  at  first, 
decreases  slightly,  and  then  increases  sharply  and  becomes  unbounded  as  t  ->  ilg .  In  the  region  where  the 
maximum  decreases,  the  pressure  disturbance  dominates  in  the  exponent  of  (3.17).  yielding  a  region  where 
the  chemical  activity  is  suppressed,  and  so  the  temperature  evolves  according  to  a  homogeneous  heat  equa¬ 
tion;  this  parallels  that  of  the  linear  pressure  case.  However,  unlike  the  linear  pressure  case,  the  pulse  does 
not  grow  unbounded,  but  instead  reaches  a  minimum  and  then  increases  back  to  zero.  Thus,  there  exists 
some  finite  time  at  which  the  pressure  no  longer  dominates  the  exponent,  ;ind  the  problem  begins  to  once 
again  evolve  according  to  classical  thermal  explosion  theory.  A  comparison  of  the  pressure  pulse  and 
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maximum  temperature  perturbation  is  provided  in  Figure  6  for  one  selected  value  of  to.  Note  that  the  max¬ 
imum  in  T j  decreases  over  most  of  the  pressure  pulse,  and  that  the  temperature  again  rises  once  the  pressure 
increases  to  some  critical  value.  Finally,  we  note  here  that  the  above  behaviour  was  seen  to  occur  numeri¬ 
cally  for  all  pressure  amplitudes  less  than  about  -2.8.  and  also  occurs  for  non-zero  values  of  ji7 . 

3.1.3.  SINUSOIDAL  PRESSURE  VARIATION  Let  the  flow  be  subject  to  a  sinusoidal  pressure 
variation  at  i  =  0  of  the  form 

Pc(t)  -  PA  sin  (or  +  d>).  (3.21) 

where  to  >  0  is  the  frequency.  PA  is  the  amplitude,  and  d>  is  the  phase  of  the  pressure  variation.  This  profile 
was  chosen  since  it  can  be  regarded  as  a  single  Fourier  mode  of  a  more  general  pressure  disturbance.  Figure 
7(a,b,c)  are  plots  of  tlg  versus  the  phase  angle  O  for  PA  -  1,  2,  4,  respectively,  and  for  various  values  of  co 
with  pT  =  0.  The  case  of  PA  =0  corresponds  to  that  of  Linan  and  Crespo  (1976)  and  is  shown  in  each 
figure  as  a  dashed  line  for  reference.  We  begin  by  noting  that  there  are  several  important  features  which 
are  common  to  all  three  figures.  These  are:  (i)  the  maximum  ignition  time  at  a  fixed  amplitude  corresponds 
to  zero  frequency  with  a  phase  tingle  greater  than  rt  (corresponding  to  an  initial  negative  pressure  wave), 
and  that  this  maximum  decreases  and  shifts  to  lower  phase  angles  as  co  is  increased  from  zero;  (ii)  for  fre¬ 
quencies  0)  =  1  and  greater,  the  ignition  time  always  lies  below  that  of  the  dashed  line  for  any  phase  angle 
and  amplitude,  indicating  that  ignition  is  enhanced  over  the  zero  pressure  disturbance  case;  and  (iii)  for  the 
high  frequency  range  (w  >  5),  the  ignition  time  is  approximately  a  constant  over  the  entire  range  of  phase 
angles,  and  that  this  constant  value  decreases  as  the  pressure  amplitude  PA  increases. 

The  same  interesting  feature  that  occurred  for  the  pressure  pulse  case  is  also  seen  to  occur  for  the 
sinusoidal  pressure  case,  in  that  the  ignition  time  undergoes  a  sudden  jump  in  its  value  for  certain  combina¬ 
tions  of  amplitude,  frequency,  and  phase.  This  jump  phenomenon  in  the  ignition  time  is  clearly  shown  in 
Figure  7c  for  PA  =  4.  For  the  curves  w  >  0.1,  the  ignition  time  drops  slowly  as  d>  is  increased  from  zero, 
reaches  a  minimum,  and  then  begins  to  increase  slowly  as  is  increased  further  until  a  critical  value  of  d> 
is  reached,  where  the  ignition  time  then  undergoes  a  sudden  jump  in  its  value.  As  d>  increases  still  further, 
the  ignition  time  reaches  a  maximum  and  then  begins  to  decrease  back  to  its  original  value  at  =  0.  To 
better  explain  this  jump  phenomenon,  we  plot  in  Figure  8  the  maximum  of  T i  versus  time  for  PA  -  4, 
(0  —  0.1  and  <t>  =  n.  The  corresponding  pressure  wave  is  displayed  as  the  dashed  curve.  In  the  region  where 
the  maximum  temperature  perturbation  is  almost  constant,  the  pressure  disturbance  dominates  in  the 
exponent  of  (3.17)  yielding  a  region  where  the  chemical  activity  is  suppressed,  and  so  the  temperature 
evolves  according  to  a  homogeneous  heat  equation,  and  ignition  is  delayed.  As  the  pressure  rises  to  some 
critical  value,  chemical  activity  resumes  and  ignition  soon  takes  place  according  to  classical  thermal  explo¬ 
sion  theory. 

Figures  9  and  10  show  that  regions  of  chemical  inactivity  can  exist  in  the  ignition  regime  even  in  the 
high  frequency  limit.  Figure  9  shows  the  maximum  of  7",  as  a  function  of  time  for  PA  =  2,  to  =  5,  and 
d>  =  n!  10,  while  Figure  10  is  a  similar  graph  except  now  to  =  10.  In  each  graph,  the  corresponding  pres¬ 
sure  variation  is  shown  as  a  dashed  curve.  Note  that  there  are  now  several  regions  where  the  maximum 
temperature  perturbation  decreases  and  hence  chemical  activity  is  suppressed,  with  the  pressure  variation 
being  negative  in  each  of  those  regions. 
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Finally,  we  note  here  that  the  general  behaviours  described  in  Figures  7-10  for  j}T  =  0  are  seen  to 
occur  for  any  value  of  pr .  The  only  exception  being  for  pr  <  0,  the  ignition  time  is  greater  than  that  given 
for  P  =  0,  while  for  pr  >  0  it  is  less. 


3.2.  R  >  0.  When  vortical  motion  is  taken  into  account,  the  governing  reactive  convection- 
diffusion  equations  for  the  temperature  and  mass  fraction  perturbations  are 

T{  ,  +  ^-—[1  -exp(-r2/4(Sc)lT,  6  =  V2T,  +<$>F[  F{  /'+Pr<'J.  (3.22) 

r 

Ffj  +  ^  [1  -  exp(-r2/4r  Sc)]Fj  e  =  V2f',  j  =  1,2,  (3.23) 

subject  to  the  conditions  (3.10).  This  system  was  solved  numerically  for  a  range  of  R ,  Pr  =  0,  Sc  =0=1, 
and  with  an  assumed  form  for  the  pressure  variation  term  Fc  (/).  Implementation  of  boundary  conditions  is 
facilitated  if  the  system  is  recast  in  Cartesian  coordinates.  The  solution  technique  is  a  2nd-order  finite 
difference  scheme  on  a  nonuniform  mesh.  To  resolve  the  structure  in  the  core  region  of  the  field,  a  coordi¬ 
nate  stretching  is  used.  To  avoid  the  singularity  at  the  origin,  no  mesh  points  are  placed  there.  The  outer 
boundaries  are  set  at  50  or  200  in  the  y -direction,  and  20  or  50  in  the  z -direction.  Grid  resolution  studies 
which  at  least  doubled  the  computational  mesh  were  carried  out  to  ensure  that  structures  were  well  resolved. 
The  resolutions  required  ranged  from  a  642  mesh  to  a  2562  mesh  for  large  vortex  Reynolds  number.  The 
time-stepping  scheme  is  a  four-stage  Runge-Kutta  which  is  formally  2nd-order  but  has  an  extended  stability 
region  making  it  accurate  and  robust  for  moderately  stiff  problems.  All  runs  were  performed  on  a  Cray 
YMP. 


As  t  increases,  the  solution  for  Tt  becomes  unbounded  at  some  finite  time  (t,g)  and  location  (y,g ,  z,g ). 
This  characterizes  the  ignition  regime.  The  special  case  R  =  0,  in  which  two  initially  unmixed  species  are 
allowed  to  diffuse  without  the  mixing  generated  by  the  vortex,  has  been  presented  in  subsection  3.1.  For  this 
case,  ignition  takes  place  along  a  line  parallel  to  the  y  -axis  which  is  located  at  zig  =  0  for  pr  =  0  and 
resides  in  the  hotter  region  for  pr  *  0.  For  R  >  0,  previous  work  [Macaraeg,  Jackson,  and  Hussaini  (1992)] 
has  shown,  in  the  absence  of  a  pressure  disturbance,  that  ignition  occurs  at  a  point  rather  than  along  an 
entire  line.  The  ignition  point  is  located  at  the  origin  for  pT  =  0  and  for  any  value  of  R .  and  resides  in  the 
hotter  region  for  Pr  *  0.  For  pr  *  0,  the  ignition  location  spirals  clockwise  towards  the  viscous  core  center 
as  R  is  increased  from  zero,  with  ignition  taking  place  within  the  core  for  R  >  70. 

To  determine  the  combined  effects  of  the  vortex  and  a  pressure  disturbance  on  the  ignition  time,  we 
produce  here  selected  results  using  the  sinusoidal  pressure  function  (3.21)  for  Pc(t).  In  all  calculations,  we 
fix  5c  =0=1  and  Pr  =  0,  and  vary  the  parameters  PA ,  R ,  to,  and  0.  Since  pr  =  0  and  R  >  0,  ignition 
will  occur  in  the  vortex  center.  Figure  ll(a,b)  are  plots  of  tig  versus  the  phase  angle  0  for  PA  =2  and  4, 
respectively,  and  for  vortex  Reynolds  number  R  =  0  (solid)  and  R  =  100  (dash)  with  to  =  0.5.  The  results 
with  R  =  0  (solid)  are  taken  from  Figure  7  and  given  here  as  a  reference  when  R  >  0.  From  Figure  1 1  we 
see  that  the  effect  of  increasing  R  decreases  slightly  the  ignition  time  for  any  phase.  This  is  consistent  with 
previous  results  obtained  when  pressure  disturbances  was  ignored  [Macaraeg,  Jackson,  and  Hussaini  (1992)]. 
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To  illustrate  the  structure  in  the  high  frequency  limit,  we  plot  in  Figures  12  and  13  the  maximum  of 
the  temperature  perturbation  7,  as  a  function  of  time  for  or  =  5  ttnd  10.  respectively,  and  for  various  values 
of  the  vortex  Reynolds  number  R  with  PA  =  2,  =  re/  10.  and  =  0.  The  special  case  of  R  =  0  was  pre¬ 

viously  shown  in  Figures  9  and  10.  respectively,  and  is  used  here  as  a  reference  case  for  when  R  >  0.  In 
both  figures,  the  structure  is  similar  to  the  R  -  0  case,  except  that  increasing  R  decreases  the  ignition  time 
slightly  until  about  R  =  100,  when  the  ignition  time  remains  constant  as  R  is  increased  further.  As  stated 
earlier,  the  temperature  maximum  always  occurs  in  the  vortex  center. 


3.3.  APPROXIMATION  FOR  THE  IGNITION  TIME,  Since  a  certain  amount  of  mixing  must 
take  place  before  ignition,  a  surprisingly  good  approximation  to  the  ignition  time  in  the  presence  of  a  pres¬ 
sure  disturbance  can  be  found  by  considering  the  "homogeneous"  problem 


dT\  1  r,  +  p  in 

- =  —  e 

dt  ,9 


(3.24) 


where  the  factor  t°  is  the  ignition  time  in  the  absence  of  a  pressure  disturbance  chosen  to  be  the  Linan  and 
Crespo  (1976)  result,  which  accounts  for  the  role  that  diffusion  plays  in  the  mixing  of  the  reactants  as  well 
as  the  role  that  the  initial  temperature  profile  plays.  For  pr  = -2,  0,  and  2,  we  find  that  if  =  13.1961. 
5.8285.  and  1.7736,  respectively,  for  our  choice  of  the  time  scales  (3.4).  Integrating  (3.24),  we  find  that 


T I 


-M 


P,U) 


he  0 


dt 


J 


(3.25) 


Ignition  takes  place  when  the  quantity  in  the  brackets  of  (3.25)  vanishes:  that  is,  the  "homogeneous"  ignition 
time  t,"  is  given  implicitly  by  the  equation 


A(t) 


d  x  =  i, 


,P 


(3.26) 


In  the  absence  of  a  pressure  disturbance,  Pc{t)  =  0  and  therefore  if,  =  if.  Below  we  give  selected  com¬ 
parisons  between  the  ignition  times  predicted  by  this  "homogeneous"  theory  and  by  the  full  numerical  solu¬ 
tions  of  the  previous  subsections. 

First,  consider  the  simple  linear  pressure  profile  given  by  (3.19).  Substitution  into  (3.26)  yields  for  the 
approximate  ignition  time 


—  In  [  1  -©/£].  (3.27) 

© 

Note  that  (3.27),  in  addition  to  defining  the  approximate  ignition  time,  also  approximates  the  critical  value 
of  ©  for  which  ignition  can  take  place.  That  is,  ignition  occurs  for  0  <  ©  <  ©/',  where  ©/'  =  1  /  if  is  the 
critical  value;  for  w  >  ©rH,  ignition  does  not  take  place.  In  Figure  14(a)  we  plot  the  ignition  times  for  the 
numerical  solutions  of  subsection  3.1.1  (solid)  and  the  approximate  solutions  (dashed)  for  (3r  =  -2.  0,  and  2 
as  a  function  of  ©.  For  each  pr ,  note  that  the  approximate  ignition  times  compare  rather  well  to  that  of  the 
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lull  solutions.  Now  consider  the  sinusoidal  pressure  function  men  by  (3.21).  For  this  case,  (3.26)  must  be 
solved  implicitly  for  t"  for  fixed  values  of  PA ,  co.  <1>  ;md  (Jr.  In  Figures  14(b.c)  we  plot  the  ignition  times 
for  the  numerical  solutions  of  subsection  3.1.3  (solid)  and  the  approximate  solutions  (dashed)  as  a  function 
of  phase  <J>  and  for  PA  =  4.  to  =  0.1  ;tnd  PA  =  4,  to  =  0.5,  respectively,  with  (i,  =  0.  In  each  figure  note  how 
weff  (he  approximate  ignition  times  compare  to  the  ignition  times  of  the  full  numerical  problem.  Thus,  we 
see  that  our  theory  of  "homogeneous"  ignition  predicts  surprisingly  well  the  true  ignition  times.  Finally,  we 
note  here  that  our  definition  of  the  "homogeneous”  ignition  time  can  also  be  used  to  predict  ignition  when 
vortical  motion  is  present,  since  ignition  times  vary  little  with  vortex  Reynolds  number  R  . 

4.  DIRECT  NUMERICAL  SIMULATION 

To  investigate  post-ignition  events  up  to  the  development  of  a  flame  sheet,  the  direct  numerical  solu¬ 
tion  to  (2.3)  is  necessary.  Unfortunately,  as  mentioned  above,  the  numerical  solution  to  the  full  system  (2.3) 
presents  a  formidable  task.  One  approximation  which  will  simplify  (2.3)  to  a  more  tractable  system  ami  has 
often  appeared  in  the  literature  is  the  assumption  of  small  heat  release  (e.g.,  Malkowsky  and  Sivashinsky 
(1979);  Van  Harten,  Kapila  and  Matkowsky  (1984)].  The  resulting  system  can  then  be  solved  numerically. 
We  begin  the  analysis  of  small  heat  release  by  assuming  p  <<  1.  The  temperature  and  density  must  now  be 
rescaled  according  to 


7*  -*  1  +  P (T  +  Pr (r )),  p^l+Pp. 

while  the  other  variables  remain  unchanged.  To  leading-order  in  p.  the  system  (2.3)  reduces  to 

7,  +  — ]1  -  exp(-/-/4r  ,Sc)|  =  V-  r  +  <{>A.  ‘F,/-  ->c 

r' 


(4.1) 


(4.2) 


Fj,i  +  ^[1  -  exp(-r2/4/  Si')\Fj'f)  =  W2Fj  -  <j>r‘  F ,  F.e™  +  P'  U'\  j  =  1,2,  (4.3) 

for  the  temperature  and  mass  fractions,  while  the  leading  order  solution  to  the  continuity  and  momentum 
equations  (2.3b)-(2.3e)  is  the  incompressible  vortex  solution  (3.5)-(3.7).  The  density  can  be  found  from  the 
gas  law.  In  the  above  equations,  X  =  pZe  is  the  product  of  the  heat  release  parameter  and  Zeldovich 
number,  8  =  Py/(Y-  1),  and  (3.4)  with  A  =  <p  has  been  used  in  the  delinition  of  the  Damkohler  number. 
Finally,  we  note  here  that  by  using  the  definition  (4.1)  for  the  temperature,  we  have  subtracted  out  the  time 
varying  part  due  to  pressure  fluctuations;  that  is,  if  T  oscillates  at  tiny  point  in  the  flow  field,  it  will  be  due 
solely  to  the  nonlinear  source  term  in  (4.2). 

The  system  (4.2)-(4.3)  must  be  solved  numerically  subject  to  the  boundary  and  initial  conditions 

T  =  0,  F\  =  1,  F2  =  0  at  i  =  0.  r  >  0,  0  <  0  <  n.  and  t  >  0,  /  -»  0  <  0  <  n.  (4.4a) 

T  =  Pr,  F |  =  0,  F2  =  <p_l  at  f  =  0.  /■  >  0.  n  <  0  <  2rt  ,  ;md  /  >  0,  r  — »  n  <  0  <  2 n.  (4.4b) 

This  system  was  solved  by  the  finite  difference  scheme  mentioned  in  Section  3.  To  illustrate  the  numerical 
solution  of  the  continuous  evolution  from  nearly  frozen  flow  to  near  equilibrium  flow,  we  produce  here 
selected  results  using  the  sinusoidal  pressure  function  (3.21)  for  F,.(t).  In  all  calculations  we  take 
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Sc  =  <}>  =  1  and  Pr  =  0. 


To  investigate  the  combined  effects  of  a  pressure  disturbance  and  of  vortical  motion  on  the  structure 
of  the  full  solution,  we  begin  by  examining  the  structure  in  the  absence  of  these  effects.  Figure  15(a)  is  a 
plot  of  the  maximum  values  of  7  as  a  function  of  time  for  R  =  0,  Pc(t)  =  0  and  X  =  1,  0.1  and  0.01.  Since 
pr  =  0  and  R  =  0,  this  maximum  occurs  along  the  line  z  =  0.  For  X  =  l,  the  temperature  increases  gradu¬ 
ally  from  zero  at  t  =  0  to  7  =  0.5  as  t  -» this  maximum  value  corresponding  to  the  flame  sheet  value 
given  in  the  next  section.  As  X  decreases  from  one,  the  rise  in  temperature  near  /  =  0  increases  sharply, 
indicating  the  existence  of  a  well-defined  ignition  regime.  An  interesting  trend  develops  if  pressure  fluctua¬ 
tions  are  now  included.  In  Figure  15(b)  we  plot  the  maximum  values  of  the  temperature  as  a  function  of 
time  for  R  =0  and  for  the  sinusoidal  pressure  function  (3.21)  with  PA  =  4,  to  =  .5  and  d>  =  1.25  n.  For 
X  =  1  and  0.1,  the  temperature  rises  from  t  =  0  to  a  maximum  value  slightly  below  the  flame  sheet  value  of 
0.5,  and  then  begins  to  oscillate  in  time  with  decreasing  amplitude;  as  /  -» °°,  the  oscillations  die  out  and 
the  temperature  maximum  approaches  the  value  0.5.  As  X  is  decreased  from  0.1  to  0.01,  ignition  takes 
place  so  rapidly  that  the  acoustic  field  is  seen  to  have  almost  no  effect  on  the  temperature  field.  The 
absence  of  any  effect  due  to  pressure  disturbances  is  consistent  with  the  flame  sheet  solution  presented  in 
the  following  section.  A  surprising  phenomenon  occurs  when  vortical  motion  is  included.  In  Figure  15(c) 
we  again  plot  the  maximum  values  of  the  temperature  as  a  function  of  time,  again  for  X  =  1,  0.1  and  0.01, 
with  R  =  50  and  the  same  pressure  profile  as  in  Figure  15(b).  Since  p7  =  0  and  R  >  0,  this  maximum 
occurs  in  the  vortex  center.  Note  that  for  all  three  values  of  X,  the  temperature  profile  evolves  in  time  with 
little  influence  from  the  pressure  field.  That  is,  the  temperature  within  the  viscous  core,  once  it  reaches  a 
maximum,  remains  almost  constant  with  time,  indicating  that  viscous  effects  damp  out  the  acoustics.  In  con¬ 
trast,  we  show  in  Figure  16  the  maximum  temperature  at  the  origin  (solid)  and  the  maximum  temperature 
far  away  from  the  core  (dashed)  as  a  function  of  time  for  the  same  conditions  as  in  Figure  15(c)  with  X  =  1. 
From  this  figure  we  see  that,  although  the  temperature  within  the  core  remains  almost  constant  with  time, 
the  outer  flow  is  oscillating.  This  is  consistent  with  Figure  15(b),  since  the  outer  flow  does  not  feel  the 
influence  of  the  vortex  so  that  the  acoustic  field  has  a  noticeable  effect.  Finally,  the  oscillations  in  the  outer 
flow  will  eventually  die  out,  consistent  with  the  flame  sheet  solution  given  in  the  next  section. 

5.  DIFFUSION  FLAME  REGIME 

After  ignition  has  taken  place,  a  thin  diffusion  flame  exists  and  is  characterized  by  a  chemical  reaction 
time  that  is  much  smaller  than  a  characteristic  diffusion  time.  Chemical  reactions  then  occur  in  a  narrow 
zone  between  the  fuel  and  the  oxidizer,  where  the  concentrations  of  both  reactants  are  very  small. 
Mathematically,  the  assumption  of  very  fast  chemical  reaction  rates  leads  to  the  limit  of  infinite  Damkohler 
number  which  reduces  the  diffusion  flame  to  a  flame  sheet  (i.e.,  local  chemical  equilibrium).  This  assump¬ 
tion  significantly  reduces  the  complexity  of  the  problem  since  it  eliminates  the  analysis  associated  with  the 
chemical  kinetics.  For  many  flows,  the  assumption  of  local  chemical  equilibrium  adequately  predicts  the 
location  and  the  shape  of  the  diffusion  flame  [Buckm aster  and  Ludford  (1982);  Williams  (1985)].  For  finite 
values  of  the  Damkohler  number,  equations  (4.2)  must  be  solved  numerically,  and  this  was  done  in  the  pre¬ 
vious  section. 

We  begin  the  description  of  the  diffusion  flame  regime  by  defining  the  following  conserved  variables 
[see,  e.g.,  Williams  (1985)]: 
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r+F,  =  pT+(1-Pr)Z. 


r  +  F2  =  (Pr  +<t>'1)(l  -Z). 


-i 


z  = 


F  i  ~  F  2  +  0' 
1+0-'  * 


where  Z  is  the  mixture  mass  fraction  and  satisfies  the  convection-diffusion  equation 

Z ,  +  ^-[1  -exp(-/-2/4/Sc)]Ze  =  V2Z, 
r 

subject  to  the  initial  and  boundary  conditions 

Z  =  1  at  t  -  0,  r  >  0,  0  <  0  <  it,  and  f  >  0,  r  ->  «,  0  <  0  <  it, 

Z  =  0  at  t  =  0,  r  >  0,  it  <  0  <  2n,  and  f>0,  r  -+<*>.  jc<Q<  2tt. 

In  the  limit  X  -+  0  (i.e.,  infinite  Damkohler  number)  the  flame  sheet  solution  is  given  by 

F,  =  l-(l  +  <|r,)(l-Z).  F2  =  0, 

r  =  (pr  +0_1)(1-Z), 


valid  for  Z  >Zf,  and 


F,  =  0,  F2  =  0-*  -  (1  +  r‘)Z, 


T  =  pr  +  (1  -  Pr)Z, 


(5.1) 

(5.2) 

(5.3) 


(5.4) 


(5.5a) 

(5.5b) 


(5.6a) 

(5.6b) 


(5.7a) 


(5.7b) 


valid  for  Z  <Zf.  Here,  Z{  defines  the  location  of  the  flame  sheet  where  both  the  reactants  vanish,  given  by 
the  implicit  relation 


and  T  takes  the  adiabatic  flame  value 


T  _  1  + 

T/  1+0 


(5.8a) 


(5.8b) 


Note  that  the  flame  location  is  independent  of  pr  and  p.  Once  Z  =  Z(r,  0,0  is  known,  then  the  other  vari¬ 
ables  ( r,F(,F2 )  can  be  found  from  (5.6)-(5.7).  We  remark  here  that  if  pr  =  0  and  0  =  1,  then  Tf  =  0.5, 
which  is  consistent  with  Figures  15-16.  We  also  note  that  the  flame  sheet  structure  and  position  is 
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independent  of  the  pressure  disturbance  Pc(t),  again  consistent  with  Figures  15-16  in  the  limit  /  -»  <». 

As  mentioned  in  the  Introduction,  equation  (5.4)  for  Z  is  exactly  Marble’s  problem  (Marble  (19X5)) 
when  ignition  events  are  ignored  (i.e.,  the  vortex  is  turned  on  at  t  =  0  with  the  initial  configuration  given  by 
the  flame  sheet),  though  derived  for  order  one  heat  release  and  with  the  constant  density  approximation,  and 
was  solved  numerically  by  Laverdant  and  Candel  (198X)  and  Rehm,  Baum,  Lozier  and  Aronson  (19X9). 

6.  CONCLUSIONS 

The  problem  of  a  flame  interacting  with  a  vortex  is  supposed  to  model  certain  fundamental  mechan¬ 
isms  of  turbulent  diffusion  flames  which  form  the  basis  of  many  propulsion  devices.  Acoustic  waves  are 
always  present  in  these  combustion  systems,  so  that  it  is  important  to  study  the  role  of  acoustics  in  such 
flame/vortex  interactions.  Having  said  this,  we  refrain  from  any  further  comment  on  the  implications  of  the 
present  study  to  practical  applications. 

The  present  study  assumes  a  one-step,  irreversible  Arrhenius  reaction  between  initially  unmixed 
species  occupying  adjacent  half-planes  which  are  allowed  to  mix  and  react  by  convection  and  diffusion  in 
the  presence  or  absence  of  a  vortex,  and  investigates  the  influence  of  the  acoustic  field  (spatially  uniform 
but  time  dependent  pressure  waves  of  small  amplitude)  on  the  ignition  time  and  flame  structure  which  are 
functions  of  vortex  Reynolds  number,  initial  temperature  difference  of  the  reactants,  equivalence  ratio  and 
Schmidt  numbers. 

In  the  case  of  zero  vortex  Reynolds  number  and  near  equal  initial  temperatures  >/  the  reactants,  this 
study  concludes  that  the  low  frequency  pressure  waves  accelerate  or  decelerate  ignition  depending  on  the 
phase;  the  pressure  waves  with  moderate  to  high  frequencies  always  enhance  ignition  and  are  essentially 
independent  of  phase.  Another  key  result  is  that  the  simple  theory  of  homogeneous  ignition  developed  in 
this  study  is  found  to  predict  these  trends  very  well,  and  therefore  it  is  proposed  as  a  good  engineering 
approximation  for  obtaining  the  ignition  times  for  the  more  complicated  flow  fields.  When  vortical  motion 
was  included,  the  overall  ignition  time  decreases  slightly  and  then  asymptotes  to  a  finite  value  with  increas¬ 
ing  vortex  Reynolds  number,  consistent  with  previous  results  obtained  by  ignoring  acoustics  [Macaraeg. 
Jackson,  and  Hussaini  (1992)].  While  the  asymptotic  analysis  is  confined  to  incipient  ignition,  a  direct 
numerical  simulation  is  carried  out  to  investigate  post  ignition  events.  This  is  limited  by  the  small  heat 
release  assumption.  Relaxing  this  rather  drastically  restrictive  assumption  is  worth  a  serious  attempt.  An 
important  conclusion  from  the  numerical  investigation  is  that  acoustics  has  little  effect  on  the  flow  field  in 
the  vortex  core,  while  the  flow  outside  the  core  oscillates.  It  is  believed  that  viscosity  in  the  reacted  core 
damps  out  the  acoustic  field. 
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Figure  5.  Plot  of  the  maximum  of  the  temperature  perturbation  Tj  as  a  function  of  time  t  for  several  values 
of  to  and  for  PA  -  -4  and  PT  =  0.  Pulse  Pressure  Profile. 
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Figure  6.  Plot  of  the  maximum  of  die  temperature  perturbation  T,  as  a  function  of  time  for  to  =  0.005. 
PA  -  -  4,  and  pr  =  0.  The  dashed  curve  is  the  pressure  profile.  Pulse  Pressure  Profile. 


Figure  7.  Plot  of  the  ignition  time  tig  as  a  function  of  phase  O  for  several  values  of  the  frequency  o)  and  for 
Pt  =  0;  7(a):  PA  -  1;  7(b):  PA  =  2;  7(c):  PA  =  4.  The  dashed  curve  corresponds  to  the  results  of  Linan 
and  Crespo  (1976/.  Sinusoidal  Pressure  Profile. 
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Figure  9.  Plot  of  the  maximum  of  the  temperature  perturbation  7",  as  a  function  of  time  for  PA  =  2.  to  =  5, 
<b  =  it  / 10,  Pr  =  0.  The  dashed  curve  is  the  corresponding  sinusoid.  1  pressure  profile. 


Figure  10.  Plot  of  the  maximum  of  the  temperature  perturbation  T\  as  a  function  of  time  for  PA  =  2, 
to  =  10,  -  n/ 10,  pr  =  0.  The  dashed  curve  is  the  corresponding  sinusoidal  pressure  profile. 
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Figure  14.  Plot  of  the  ignition  times  corresponding  to  the  full  numerical  solutions  (solid)  and  the  homo¬ 
geneous  ignition  times  (dashed)  as  a  function  of  (a)  a)  using  the  linear  pressure  profile:  (b)  phase  d>  using 
the  sinusoidal  pressure  profile  with  PA  =  4  and  to  =  0.1;  and  (c)  phase  d>  using  the  sinusoidal  pressure 
profile  with  PA  =  4  and  ©  =  0.5. 
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